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Introduction 

The howardites, eucrites, and diogenites (HEDs) 
are a suite of related meteorite types that formed by 
igneous and impact processes on the same parent body 
[1], Multiple lines of evidence, including infrared 
spectroscopy of the asteroid belt and the petrology and 
geochemistry of the HEDs, suggest that the asteroid 4 
Vesta is the parent body for the HEDs [e.g., 2, 3]. Ob- 
servations by NASA’s Dawn spacecraft mission 
strongly support the conclusion that the HEDs are 
from Vesta [e.g., 4, 5]. 

The abundances of the moderately siderophile ele- 
ments Ni, Co, Mo, W, and P in eucrites require that 
most or all of the metallic phase in Vesta segregated to 
form a core prior to eucrite solidification [6-8]. These 
observations place important constraints on the mode 
and timescale of core formation on Vesta. Possible 
core formation mechanisms include porous flow, 
which potentially could occur prior to initiation of sili- 
cate melting, and metallic rain in a largely molten sili- 
cate magma ocean. Once the core forms, convection 
within the core could possible sustain a magnetic dy- 
namo for a period of time. We consider each process in 
turn. 

Core Formation: Porous Flow 

Melting in the Fe-FeS system begins at the eutectic 
temperature of 988 °C at 30 weight % S. Inclusion of a 
small amount of Ni will lower the eutectic temperature 
slightly and shift it towards lower S [9]. Existing mod- 
els of the thermal evolution of Vesta all assume an Fe- 
Ni-S solidus of 940 °C and a near-eutectic liquidus of 
960 °C [10-12] and thus require a relatively large 
amount of S in Vesta’s core. 

Assuming Vesta has sufficient S to have near- 
eutectic melting of the metal phase, the metal would be 
completely molten prior to the onset of silicate melt- 
ing. At low melt fractions, liquid iron in solid silicates 
has a high dihedral angle and forms isolated pockets of 
melt [13]. The dihedral angle is a function of oxygen 
fugacity, with high values of f0 2 permitting connected 
sheets of liquid Fe [14], but Vesta likely formed under 
low oxygen fugacity conditions (at or below the iron- 
wustite buffer, [e.g., 7]). Even with a high dihedral 
angle, once the melt fraction of metal exceeds a few 
percent, connected melt channels are possible [15, 16]. 
Because of the low viscosity of liquid iron metal and 
its high density relative to silicates, Darcy flow in such 
connected melt channels will be efficient, allowing 
geologically rapid separation of the liquid metal from 


the solid silicate. However, this separation is likely to 
be imperfect: once most of the metal has drained from 
the silicate and the melt fraction drops below the criti- 
cal level for maintaining connected channels, the chan- 
nels may pinch off, trapping the last few percent of 
metal in the solid. 

Core Formation: Magma Ocean Metallic Rain 

A significant problem with the near-eutectic melt- 
ing porous flow model outlined above is that the eutec- 
tic composition requires a large amount of sulfur (~30 
weight %, [9]). In contrast, a recent study synthesized 
a broad range of constraints on Vesta’s composition 
and favored a model composed of 75% H chondrite 
and 25% CM chondrite [17]. This corresponds to a 
core of 16 % S, 73.7 % Fe, and 10.3 % Ni, which is far 
removed from the Fe-FeS eutectic composition. For 
this low S abundance, melting of metal still begins at 
the eutectic temperature but does not conclude until 
about 1350 °C [9]. Because the HEDs are relatively 
depleted in volatile elements, it is possible that the S 
abundance will be even lower than this, which would 
shift the Fe-FeS liquidus to an even higher tempera- 
ture. The solidus for the primitive silicate H+CM man- 
tle is likely in the range 1100-1150 °C [18, 19]. Thus, 
metal and silicate melting are likely to overlap in time 
during the formation and differentiation of Vesta. It is 
worth noting that all published Vesta-specific thermal 
models assume that metallic melting is complete well 
before the onset of silicate melting [10-12]. Although 
the overall picture derived in those models is probably 
generally correct, some of the details of the inferred 
thermal evolution and differentiation history in those 
models must be incorrect. For example, in order to use 
Hf-W radiometric ages for metal-silicate separation 

[20] as a constraint on the thermal models, the metal 
liquidus used in the thermal model must be correct. 

The viscosity of a suspension of silicate crystals in 
a silicate melt decreases by many orders of magnitude 
when the liquid fraction increases from 40% to 60% 

[21] . When full melting of the Fe-Ni-S metal phase is 
achieved at 1350-1400 °C, the silicates will be near 
50% melt fraction. Assuming that Vesta formed in the 
first 2-3 half-lives of 26 A1 after formation of CAIs in 
the solar nebula and incorporating the insulating ef- 
fects of a low conductivity megaregolith layer [22], 
radioactive decay of 26 A1 can easily heat Vesta to > 
1400 °C, including the effects of latent heat for both 
the silicate and metal phases. At this stage of evolu- 
tion, Vesta would effectively be a magma ocean. 


Given the low viscosity of the silicate melt and the 
high density contrast between the metal and the sili- 
cate, metal blobs will fall to the center of Vesta by 
Stokes flow in a sort of metallic rain [23, 24]. Based 
on the various chemical and physical constraints out- 
lined here, we favor a magma ocean and iron metal 
rain as the most likely method of forming a core on 
Vesta. Solidification of a magma ocean can also ex- 
plain the compositional characteristics of the eucrites 
and diogenites [25]. 

Core Convection and a Magnetic Dynamo 

Definitive observations of an intrinsic magnetic 
field on Vesta would constrain the structure and rate of 
cooling of the core and thus would provide important 
constraints on Vesta’s thermal evolution. Unfortu- 
nately, Dawn did not include a magnetometer as part 
of its science payload. However, a remnanent magnetic 
field of at least 2 microtesla has been measured in the 
eucrite Allan Hills 81001 [26], which in turn implies 
that Vesta at some point in its history had a convecting 
liquid metal core. 

The short-lived radioactive isotope 26 A1 is a potent 
source of radioactive heating in the first few million 
years of Vesta’s history and is concentrated in the sili- 
cate portion of Vesta. Thermal evolution models show 
that 26 A1 heating of Vesta’s mantle initially acts as a 
thermal blanket and strongly suppresses heat flux out 
of the core. Cooling of Vesta after 26 A1 becomes ex- 
tinct eventually allows the core to briefly convect, but 
in that model core convection ends before the crust 
cools below the Curie temperature, such that evidence 
for a magnetic dynamo would not be preserved [27]. 

However, several factors may modify this conclu- 
sion. First, melt migration in the silicate magma ocean 
transfers latent heat of melting from the interior to the 
near surface. Second, magmatism also transports alu- 
minum in plagioclase from the mantle to the crust [28]. 
Both effects act to cool the mantle relative to models 
that do not include magmatism and thus reduce the 
thermal blanketing of the core. Third, 60 Fe is a possible 
radioactive heat source with a half-life that is several 
times longer than 26 A1. The abundance of 60 Fe in the 
solar nebula is controversial and was possibly spatially 
heterogeneous [29, 30], but 60 Ni anomalies in the 
eucrites Bouvante and Juvinas [31] point to the pres- 
ence of 60 Fe in Vesta. Prior thermal evolution studies 
have often either neglected 60 Fe altogether [27] or as- 
sumed that it is 100% partitioned into the core [10, 11, 
32]. It is likely that roughly equal amounts of Fe will 
be in the metal and the silicate portions of Vesta [17], 
but the much greater concentration of Fe in the core 
will contribute to driving core convection. We are cur- 
rently assessing the effects of these factors using a 
modification of our model for the thermal evolution 


and dynamo activity on Mars [33]. This model ac- 
counts for the thermal evolution of both the core and 
mantle, including magmatic heat transport, and calcu- 
lates transport of heat producing elements from the 
mantle to the crust using appropriate partition coeffi- 
cients. Compositional convection in the core may also 
play a role in driving a geodynamo [34], although top- 
down solidification of the core could create a stably 
stratified core without a dynamo [35]. 
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